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Abstract 

We present a new calculation of long-distance contributions to the K — > tttt am- 
plitudes using the 1/Nc expansion within the framework of chiral perturbation theory. 

Along these lines we compute the chiral loop corrections to the operators Qi, Q2, Qe, 

(1/2) (3/2) 

^ ' and Qs- In the numerical analysis we present values for J5g ' ' and ' ' . We discuss 

■ also the imphcations for a theoretical understanding of the AI =1/2 selection rule. 

(N 

\0 , 

O '. 1 Introduction 
00 
a^ 

Within the standard model the calculation of the K ^ mi decay amplitudes is based 
p 1' on the effective low-energy hamiltonian for AS* = 1 transitions |l|, 

^ i=l 

^ : Ciifx) = Ziifx) + ryiifx) , r = -^t/^u , = V*,Vga , (2) 

where the Wilson coefficient functions Cj(/i) of the local four-fermion operators Qi{n) 
are obtained using perturbative renormalization group techniques 0, 0. Q1-Q5 and Qj 
are composed of products of two currents, Qq and Qs are density-density operators. In 
view of phenomenological applications Qi and Q2, as well as, Qq and Qs are of special 
importance, where 

Ql = ^SLl^dL UlI^.Ul , Q2 = "^SlI^Ul UL'Jf.dL , (3) 

Qe = -8 ^ SLQugRdL, Qs = -12 ^ CgSLgRgndL, (4) 

q=u,d,s q=u,d,s 



*Talk presented at the XVI Autumn School and Workshop on Fcrmion Masses, Mixing and CP 
Violation, Lisboa, Portugal, 6-15 October 1997. 



with ql ji = |(1 ± 75)g and Cg = (2/3, —1/3, —1/3). The hadronic matrix elements 
of the operators, {Qi{fi))i = {tttt, I\Qi{fi)\K^), are difficult to calculate but can be 
estimated using non-perturbative techniques generally for /i around a scale of 1 GeV. 
In the following analysis we will use the 1/Nc expansion within the framework of the 
chiral effective lagrangian for pseudoscalar mesons. 

The calculation of chiral loop effects motivated by the l/N^ expansion was intro- 
duced in Ref. @] to investigate the A/ = 1/2 selection rule. Experimentally one finds 
Re Ao/Re A2 — 22.2 for the ratio of the isospin amplitudes. Keeping only the dominant 
terms the theoretical expression for this ratio reads 

Re Aq ^ ^i(/i)(Qi(/i))o + 2:2(/i)(Q2(/i))o + ze{n){QG{n))o 
ReA2~ zi(/i)(Qi(/i))2 + 2;2(/u)(Q2(yu))2 

It was recognized, long ago, that the vacuum saturation approximation (VSA) 
completely fails in explaining the A/ = 1/2 rule. The authors of Ref. Q considered 
chiral loop corrections to the current- current operators, Qi and Q2, and included the 
gluon penguin operator Qq at the tree level. As a net result they obtained a large 
enhancement of the A/ = 1/2 decay amplitude which, however, was found to be still 
below the experimental value {ReAo = 3.33 x 10"'' GeV). One might note that the 
agreement with experiment is not improved by including the NLO values for the Zi [0. 

The A/ =1/2 selection rule is governed by Qi and Q2, while the ratio e'/e which 
measures the direct CP violation in K titi decays is dominated by the density-density 
operators Qq and Qs, that is to say, approximately^ 

e'/e oc ?/6(Q6)o - 1/8 (Qs) 2 + 'vr- 77-77' mixing' . (6) 

It is of special importance to investigate whether the 1/A^c (chiral loop) corrections 
significantly affect the large cancellation between the gluon and the electroweak penguin 
contributions obtained at the tree level in Ref. ||^. 

In this talk we present a new analysis of the hadronic matrix elements of the op- 
erators Qi, Q2, Qq, and Qs in which we include the ffist order corrections in the 
twofold expansion in powers of external momenta, p, and the ratio 1/Nc. Our main 
improvement concerns the matching of short- and long-distance contributions to the 
decay amplitude, by adopting a modified identification of virtual momenta in the chiral 
loops. In Section 2 we briefly review the framework of the effective low-energy approach 
and discuss the matching of short- and long-distance contributions. In Section 3 we 
list the analytic expressions for the hadronic matrix elements. Finally, in Sections 4 
and 5 we present our numerical results. 

We will focus here on the computation of the K ^ irn matrix elements. A more 
detailed discussion of the general method and the operator evolution at long distances 
is presented by T. Hambye and P. Soldan in these proceedings, and some additional 
details can be found in Refs. |P, [ToU . 



^Recent reviews on CP violation in K ^ tttt decays can be found in Ref. 0. 
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2 General Framework 



We start from the chiral effective lagrangian for pseudoscalar mesons which involves an 
expansion in momenta where terms up to O(p^) are included |[Tl[|. Keeping only terms 
which contribute to the K ^ -k-k matrix elements and are leading in it reads 



+rL5{D^U^D^'U{M^U + U^M)) + rLs{M^UM^U + MU^MU^) , (7) 



with Df,U 



ir^ + iUl^, (A) denoting the trace of A and M = diag{mu, rrifi, mg). 
f and r are free parameters related to the pion decay constant and to the quark 
condensate, respectively, with r = — 2(gg)//^. The complex matrix f/ is a non-linear 
representation of the pseudoscalar meson nonet: 



U 



exp -n 
/ 



In terms of the physical states 



n 



(AaAfe) = 25, 



ab 



(8) 



h7]' 



V2K+ 



\arl 



(9) 



with 



cos6' - y2sin6' , v^fe = sin 6* + ^2 cos 6^ , 



(10) 



and the rj — rj' mixing angle 6 for which we take 6 = —19° [|T2|. a is the instanton 
parameter which reflects the explicit breaking of the Ua{^)- The symmetry violating 
term must be included when treating the singlet tjq as a dynamical degree of freedom. 

The bosonic representations of the quark currents and densities are defined in terms 
of (functional) derivatives of the chiral action and the lagrangian, respectively: 



6S 



:iii 



5 Mi 



~r(^U^ + L^df^U^d^UU^ + 2rLsU^MU^ + rH2M^) , (12) 

V 4 / ji 



and the right-handed expressions are obtained by parity transformation. Eqs. (|TT]) and 
([T^ ) allow us to express the operators Qi in terms of the pseudoscalar meson fields. 

The 1/Nc corrections to the matrix elements {Qi)i are calculated by chiral loop 
diagrams. The diagrams are ultraviolet divergent and are regularized by a finite cutoff. 
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This procedure, which was introduced in Ref. 0], is necessary in order to restrict the 
chiral lagrangian to the low-energy domain. Since we truncate the effective theory to 
pseudoscalar mesons, the cutoff has to be taken at or, preferably, below the C(l GeV). 

To calculate the amplitudes we follow the lines of Ref. 0] and identify the ultraviolet 
cutoff of the long-distance terms with the short- distance renormalization scale fi. In 
carrying out this matching we pay special attention to the definition of the momenta 
inside the loop. This question must be addressed because the loop integrals, within 
the cutoff regularization, are not momentum translation invariant. 

In the existing studies of the hadronic matrix elements the color singlet boson con- 
necting the two densities or currents was integrated out from the beginning H, [13|, |14| . 
Thus the integration variable was taken to be the momentum of the meson in the loop, 
and the cutoff was the upper limit of its momentum. As there is no corresponding 
quantity in the short-distance part, in this treatment of the integrals there is no clear 
matching with QCD. The ambiguity is removed, for non-factorizable diagrams, by con- 
sidering the two densities to be connected to each other through the exchange of the 
color singlet boson, as was already discussed in Refs. |jlO|, |15|, |16], 0, [1^. A consis- 
tent matching is then obtained by assigning the same momentum to the color singlet 
boson at long and short distances and by identifying this momentum with the loop 
integration variable (see Fig. 1). The effect of this procedure is to modify the loop 
integrals, which introduces noticeable effects in the final results. More important it 
provides an unambiguous matching of the 1/Nc expansion in terms of mesons to the 
QCD expansion in terms of quarks and gluons. 




short-distance 



long-distance 



Fig. 1. Matching of short- and long-distance contributions. 



The approach followed here leads to an explicit classification of the diagrams into 
factorizable and non-factorizable: there is no way to establish a unique momentum 
prescription for the factorizable terms. However, factorizable loop diagrams refer to the 
strong sector of the theory and give corrections whose scale dependence is absorbed in 
the renormalization of the chiral effective lagrangian. This statement is obvious in the 
case of the (conserved) currents and was proved explicitly, within the cutoff formalism, 
in the case of the bosonized densities.^ Consequently, in the factorizable sector the 
chiral loop corrections do not induce ultraviolet divergent terms, i.e., the only remaining 
ultraviolet structure of the matrix elements is contained in the overall factor ~ 
whose scale dependence exactly cancels the factorizable evolution of the density- 

^For details see Ref. Q and the contributions of T. Hambye and P. Soldan to these proceedings. 



4 



density operators Qq and Qs at short distances. The finite terms arising from the 1/Nc 
corrections to these operators, on the other hand, will not be absorbed completely in 
the renormalization of the various parameters. Nevertheless, as the complete sum of the 
factorizable diagrams is finite, we are allowed to send the scale of the factorizable terms 
to arbitrary high values, i.e., we are allowed to calculate the remaining finite corrections 
in dimensional regularization, which is momentum translation invariant and yields an 
unambiguous result for the factorizable matrix elements. Finally, the non-factorizable 
loop diagrams have to be matched to the Wilson coefficients and should cancel scale 
dependences which arise from the (non-factorizable) short-distance expansion. 



3 Hadronic Matrix Elements 



Expanding the lagrangian of Eq. (^, as well as, the densities and currents of Eqs. (|TT]) 
and (|12D in terms of the pseudoscalar fields, the various factorizable and non-factoriz- 
able contributions to the K ^ Tin matrix elements can be calculated from the diagrams 
of Figs. 2 and 3, respectively. 






(ij)(kl) 




® 

(ijXkl) 



^0 \Q^r 

(ij)(ki) (ij)(ki) 



(ij)-^>(ki) 



(ij)(kl) 



Fig. 2. Factorizable diagrams for the matrix elements of the operators Qi. Crossed circles represent the 
currents or densities, black circles denote strong vertices. The lines represent the pseudoscalar mesons. 





(ij)(kl) 



+ (ij)*Xkl) 



Oj)(W) 



(ij)(U) 



Fig. 3. Non-factorizable diagrams for the matrix elements of the operators Qi 



In addition, the wave function renormalization must be included, and the values 



of the various parameters in the lagrangian must be expressed in terms of physical 
quantities. This was done using dimensional regularization as discussed in the previous 
section. We verified that the scale dependence arising from the factorizable diagrams 
is completely canceled by the scale dependence of the various parameters in the tree 
level expression of the matrix elements. Finally, we obtain the following expressions: 



l{'K + Tl-\Ql\K^) 



0\F 



t{7r'n'\Q2\K') 



0\F 



0, 



^(7r0 7r°|gi|K^) 



0\F 



-t{7r^n-\Q2\K'') 



0\F 



X 1 + 



F2 



0\F 



t{n^n-\Qe\K") 



0\F 



2m\ 



\rh + rrii^ 



0\F 



2m^ 



K 



(13) 
(14) 

(15) 
(16) 



0\F 



'-V2F^ 



16 



m + m 

2 



1 + 



3Fr 



Iml - llml)Ll 



K 



2mi)Ll 



(17) 



with X = V2F^\ 



mt- — mz 



''K 

Fk 
F. 



and the couplings Lg, Lg defined through the relations 



ra 



K 



1 + ^« 



rh + nis 



2m 



[m 



K 



m 



F2 



— {LI-2L 



(19) 



For the reason of brevity, in the analytic expressions of this section we omit scale- 
independent (finite) terms resulting from the loop corrections which, nevertheless, are 
taken into account within the numerical analysis. Note that in Eqs. ([T^)-(|T^ we 
used 1/Fj^ and 2m^/(m + m^) rather than the bare parameters 1/f and r. Formally, 
the difference represents higher order effects. Nevertheless, the appearance of / or r 
in Eqs. (|Iip-(|I^) would induce a scale dependence which has no counterpart at the 
short- distance level and will be absorbed at the next order in the chiral expansion. 

The non-fact or izable diagrams of Fig. 3 are calculated by introducing a cutoff reg- 
ulator Ac. Using the matching prescription of Fig. 1 we obtain 



t{n"n"\Q,\K") 



0\NF 



0, 



i{7r+7r- \Q,\K^) 



0\NF 



X 



t{7r"n"\Q2\K") 



0\NF 



l{TT^Tr-\Q2\K") 



0\NF 



IQtt^F^ 
X 

X 
167r2F2 



-3Al + -{ml + 12ml)\nAl 



-Al + - (m 



6m!) In A: 



K 



2 " 



InA^ 



(20) 
(21) 

(22) 

(23) 
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^{.'.'mKr' = ^{n^7.-mKT' = tA^U-^Y^^^I (24) 

loTT'^r^A \m + rris J 

^{n'n'\Qs\Kr' = 77^2 'A ^^^l (25) 

IbTT'^l'^A \m + nis J 

^(Tr+Tr-IQslK")^^ = JI^^LV «inA2 (26) 

The scale dependent terms in Eqs. (pT])-(p^) have to be matched to the Wilson 
coefficient functions. To this end we use the numerical values presented in the leading 
logarithm analysis of Ref . |p . 

4 and 

To evaluate the hadronic matrix elements of Qq and we pursued the following strat- 
egy. First, the non-factorizable contributions were calculated, in the isospin limit, from 
the diagrams of Fig. 3. In this part of the analysis the finite terms, which are system- 
atically determined by the momentum prescription of Fig. 1, are completely taken into 
account. By using algebraic relations all integrals resulting from the diagrams of Fig. 3 
can be reduced to three basic integrals which were calculated in the cutoff formalism. 
Second, the finite terms arising from the factorizable loop diagrams of Fig. 2, as well 
as, from the renormalization of the wave functions, the masses and the low-energy 
couplings were estimated using dimensional regularization, as discussed at the end of 
Section 2. 

The values we obtain for the matrix elements of Qq and Qg are presented in Ta- 
ble 1, where we also specify the various contributions coming from the factorizable 
and the non-factorizable diagrams, respectively. In the numerical analysis we used 
the ratio ms/rh = 24.4 ± 1.5 [1^. In the matrix elements we extracted the factor 



= [2m|./(m + mg)]^, whose dependence on the (factorizable) scale is canceled ex- 
actly, for a general density- density operator, by the diagonal evolution of the Wilson 
coefficients. 

In Table 2 we list the corresponding values for the Bi factors which quantify the 
deviation of the hadronic matrix elements from the VSA results: 

Bi'^'^ = |(Q6)o/(g6)r^| , 4'/'^ = |(g8)2/(Q8)r^| . (27) 

The VSA expressions for the matrix elements were taken from Eqs. (XIX. 16) and 
(XIX.24) of Ref. §]. Numerically, they are |(Q6)o| = 35.2 ■ i?^ MeV and |(g8)2| = 
56.6 MeV ■ [i?^ — (0.389 GeV)^]. The second term in the expression for Qs contributes 
at the 2 % level and has been neglected in Table 2. 

We discuss next the corrections to the matrix elements {Qg)o and {Q8)2- As already 
mentioned, the operator Qq is special because the 0{p^) tree level matrix element is 
zero due to the unitarity of the matrix U. Nevertheless the one-loop corrections to 
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Ac = 0.6 GeV 


Ac = 0.7 GeV 


Ac = 0.8 GeV 


Ac = 0.9 GeV 


• //^ \ Lice -T lUUPo 

■ ir\ \ NF loops 


-68.4 - 37.0i 
29.8 + 37.0i 


-68.4 - 37.0i 
34.6 + 37.0i 


-68.4 - 37.0i 
39.0 + 37.0i 


-68.4 - 37.0i 
42.9 + 37.0i 


1 / /~\ \ 1 t,ot?il 
|(Q6)o| 


38.6 


33.7 


29.4 


25.5 


. / \ troo -1- P loons 
■ / ^ \ NF loops 


56.0 - O.H 
-20.7- 11. 5i 


56.0 - O.lz 
-24.8 - 11.5i 


56.0 - O.H 
-28.8 - 11. 5i 


56.0 - O.lz 
-32.8- 11. 5i 




37.2 


33.2 


29.5 


25.9 



Table 1: Hadronic matrix elements of Qe and Q% (in units of B?- MeV), shown for various values of 
the cutoff A,,. 





Ac = 0.6 GeV 


Ac = 0.7 GeV 


Ac = 0.8 GeV 


Ac = 0.9 GeV 


r(1/2) 


1.10 


0.96 


0.84 


0.72 


^(3/2) 


0.66 


0.59 


0.52 


0.46 



Table 2: Be. and -Bs factors for various values of A^. 



this matrix element must be computed. These corrections are of 0{p^ /Nc) and are 
non-vanishing. We have shown in Eqs. (pTD-(p6D that the explicit calculation of the 
loops yields a cutoff dependence (i.e., a non-trivial scale dependence) from the non- 
factorizable diagrams which has to be matched to the short-distance contribution. 
We note that in the twofold expansion in and p^/N^ the contribution of the loops 
over the 0{p^) matrix element must be treated at the same level as the leading non- 
vanishing tree contribution proportional to L5. This is revealed by the large size of 
the non-factorizable 0{p^/Nc) corrections presented in Table 1. It is the sum of both, 
the factorizable and the non-factorizable contributions, which renders the numerical 
values for {Qq)o close to the VSA value. The main effect of the loop corrections is to 
change the dependence of {Qe)o on Ac, from a flat behaviour at the tree level to the 
dependence presented in Tables 1 and 2. We note that at Ac ~ 700 MeV the value for 
the matrix element of Qq is very close to the VSA result leading to i?6 — 1- 

The Qs operator is not chirally suppressed, i.e., its C(p°) tree level matrix element 
is non-zero. The one-loop corrections, even though suppressed by a factor 1/Nc with 
respect to the leading tree level, are found to be large and negative, leading to the small 
values for Bg presented in Table 2. These large corrections persist in the octet limit 
[i.e, in the absence of the rjo, with a = b = 1 and = {Am\ — m^)/3 ]. In comparison 
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with {Qq)o, the non-factorizable corrections to (Q8)2 are less pronounced, as expected 
from the power counting scheme in p"^ and However, because the factorizable 

0{p^) and 0{p^/Nc) corrections are small (and negative), the non-factorizable terms 
produce a significant reduction of {Qs)2- The size of the higher order terms indicates 
that the leading-A'^c calculation or the closely related VSA are not sufficient for the 
matrix elements of the Qs operator. 



It is interesting to compare our results with those of other calculations. Refs. 13 



and investigated 1/Nc corrections to the matrix elements of Qq and Qs- This 
calculation considered the product of the two densities without the color singlet boson 
and the matching of short- and long-distance contributions was not explicit as in the 
present analysis. The numerical results showed also a tendency of reducing {Qs)2 
substantially, whereas (Q6)o was found to be enhanced compared to the VSA result. 
Calculations in lattice QCD obtain values for Bq close to the VSA, i?g^''^''(2 GeV) = 
1.0±0.2 [H, P and 0.76(3)(5) pf. Recent values reported for are 5f ^^^(2 GeV) = 
0.81(3)(3) 0.77(4)(4) [^, and 1.03(3) [g^]. The chiral quark model ^ yields a 
range for Bq which is above the VSA value, i?g^^^^(0.8 GeV) = 1.2 — 1.9, and predicts 
a small reduction of the Bg factor, _Bg^''^'*(0.8 GeV) = 0.91 — 0.94. Although the scales 
used in the lattice calculations and the phenomeno logical approaches are different, the 
various results can be compared as the Bq and -Bg parameters were shown in QCD to 
depend only very weakly on the renormalization scale for values above 1 GeV 0. 
We note that our result for b9^^'^^ is in rough agreement with those of the various 

(3/2) 

studies quoted above, whereas the value we obtain for B^ lies below the values 
reported previously. It is desirable to investigate whether this substantial reduction 
of {Qs)2i which is due to the non-factorizable 1/Nc corrections to the leading term in 
the chiral expansion of Qg, will be affected by higher order corrections. This point is 
of great phenomenological interest because a less effective cancellation between the Qq 
and Qg operators, in the range obtained in the present analysis, will lead to a large 
value of e' /e in the ball park of ~ 10~^. This can be confirmed or disproved by the 
forthcoming experiments at CERN (NA48), Fermilab (E832) and Frascati (KLOE). 



5 A/ =1/2 Rule 

Along the lines discussed in Section 2 we computed also the hadronic matrix elements 
of the current-current operators Qi and Q21 including the finite terms omitted in the 
analytic expressions of Section 3. Subsequently, we carried out the matching with 
the short-distance coefficient functions Zj(/i). Our numerical result for the amplitude 
Aq is depicted in Fig. 4. It shows an additional enhancement compared to the result 
of Ref. which renders the amplitude in good agreement with the observed value. 
The new contribution arises from the Qi and Q2 operators. It is due to the modified 

•^Note that due to current conservation no finite terms arise from the factorizable loop diagrams of 
Fig. 2. 
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matching prescription in the non-factorizable sector [except for approximately 20 % of 
it explained by the choice of the physical value in Eqs.(pID-(p^)]. 



MeV 




0.25 - 



0.2 - 

- 

500 550 600 650 700 750 800 850 900 950 1000 

Aj. (MeV) 

Fig. 4. ReAo in units of 10~^ MeV for 771^(1 GeV) = 150 MeV as a function of the matching scale Ac. 
The experimental value is represented by the dashed line. 

Our result is widely stable with respect to the matching scale. The main uncertainty 
displayed in Fig. 4 originates from the dependence of the Wilson coefficients on Aqcd- 
On the other hand, the A/ = 3/2 amplitude depicted in Fig. 4 is highly unstable. The 
large uncertainty can be understood from the fact that it is obtained from the difference 
of two large amplitudes of approximately the same size [namely A{K^ — >■ 7r"'"7r~) and 
A{K^ — > 7r°7r°)]. Consequently, the AI = 3/2 amplitude is not well reproduced except 
that, and this is an important point, it turns out to be sufficiently suppressed whatever 
the particular chosen scale is between 500 MeV and 1 GeV. 

In conclusion, it is certainly premature to say that the dynamical mechanism behind 
the AI = 1/2 selection rule for K ^ irn decays is now completely understood since the 
1/Nc expansion we use is only an approximate realization of non-perturbative QCD. In 
particular, vector mesons and higher resonances should be included in order to make it 
possible to choose a higher value for the matching scale (it is probable that the vector 
mesons play an important role for A2). Nevertheless we believe that the enhancement 
reported here is a further important indication that the 1/Nc approach can account for 
the bulk of the AI = 1/2 rule for K ^ nn decays. 



10 



MeV 



xlO . 




0.025 - 



r 

500 550 600 650 700 750 800 850 900 950 1000 

(MeV) 

Fig. 5. Re A2 in units of 10~^ MeV for ms(l GeV) = 150 MeV as a function of the matching scale Ac. 
The experimental value is represented by the dashed line. 
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